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A space-charge model is applied to describe the equilibrium effects of segregation of double-donor oxygen va- 
cancies to grain boundaries in dry and wet acceptor-doped samples of the perovskite oxide BaZrOa. The grain 
boundary core vacancy concentrations and electrostatic potential barriers resulting from different vacancy 
segregation energies are evaluated. Density-functional calculations on vacancy segregation to the mirror- 
symmetric S3 (112) [110] tilt grain boundary are also presented. Our results indicate that oxygen vacancy 
segregation can be responsible for the low grain boundary proton conductivity in BaZrOs reported in the 
literature. 
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While the benefits of ceramic materials make solid ox- 
ide proton conductors like barium zirconate (BaZrOs, 
BZO) desirable for use as electrolytes in electrochemical 
devices such as hydrogen fuel cells, limitations in proton 
conductivity have thus far prevented successful imple- 
mentation.^'^ It has become apparent that boundaries 
between grains in the material are the prime source of 
inhibited overall proton conductivity. '^^^ 

Experiments have established that the grain boundary 
(GB) proton resistivity is an intrinsic effect, not caused 
by the segregation of secondary phases at the GB.^'^ Two 
different explanations have instead been put forward^ , in 
which the GB effects are believed to originate in either 
a structural distortion in the GB region'^, or in the ap- 
pearance of positively charged GBs, caused by a change 
in chemical composition and leading to Schottky barriers 
and the depletion of mobile protonic charge carriers. ^'^'^'^ 
At present the latter explanation model is the predomi- 



nating one' 



6,9,10 



but the details of the GB cores are nei- 



ther well understood nor sufficiently explored. 

In this letter we present an investigation concerning 
the equilibrium segregation of oxygen vacancies to GBs 
in dry and hydrated acceptor-doped BZO and the electro- 
static potential barrier such segregation causes. Density- 
functional theory (DFT) is used to evaluate the energy of 
vacancy segregation to a E3 tilt GB and a space-charge 
model is applied to describe the equilibrium barrier and 
core vacancy concentration resulting from different segre- 
gation energies. We show that segregation energies on the 
order that we find can cause potential barriers consistent 
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with experimental findings on GB proton conductivity in 
BZO. 

The system is modeled using the periodic supercell 
technique and all DFT calculations are performed within 
the plane-wave approach as implemented in the Vienna 
ab-initio simulation package (VASP)^^'^^. Electron-ion 
interactions are described by the projector augmented 
wave method"'^'^ and a generalized gradient approximation 
(GGA) exchange-correlation functional due to Perdew 
and Wang (PW91)^^ is used. More computational de- 
tails have been reported elsewhere. In the present work 
all calculations are performed non-spin-polarized with a 
plane-wave cutoff of 400 eV (constant volume) or 520 eV 
(volume relaxations), and a 4 x 2 x 1 Monkhorst-Pack 
/c-point grid is used to sample the Brillouin zone. 

Theoretical and experimental work on the currently 
considered E3 (112) [110] tilt GB has previously been 
performed on strontium titanate^^ and this particu- 
lar GB is well suited for DFT calculations. It con- 
sists of alternating BaZrO and O2 planes and here we 
consider the mirror-symmetric BaZrO terminated case, 
which has the lowest energy. The interplanar distance is 
d = ao/-\/24 = 0.87 A, with a computed value of the bulk 
lattice constant oq — 4.25 A. Our supercell has the di- 
mensions (-s/Soo, 2-\/2ao, 3-\/6ao) in the directions ([111], 
[110], [112]) and consists of 180 atoms with the distance 
l%d = 15.6 A between the periodically repeated GBs. 

To begin with, the system is optimized by relax- 
ing the atomic structure, resulting in the configuration 
shown in Fig. 1. The GB expands 0.135 A in the per- 
pendicular direction and the calculated GB energy^^ is 
CGB = 0.78 i/ir? . Next, an oxygen vacancy is intro- 
duced. Since double-donor vacancies arc dominating in 
acceptor doped BZO, we consider the q ~ +2 charge state 
only. The vacancy formation energy AE'gg is determined 
as function of vacancy position, both with and with- 
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FIG. 1. The relaxed structure of the mirror-symmetric E3 
(112) [110] tilt grain boundary. 
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FIG. 2. Segregation energies of oxygen vacancies in the S3 
(112) [110] tilt grain boundary as function of the location of 
the vacancy, calculated using DFT. 

out relaxing the structure. As reference, the vacancy 
formation energy for the bulk system AE^^^^^^ with the 
same size of the superccU is also determined. The seg- 
regation energy is then obtained by taking the difference 



A^buii<- Note that the error due to in- 



teraction between charged defect images in the periodic 
supercell calculations is to a large degree cancelled for 
the quantity i^scg- 

We find a significant segregation tendency to layer 1 of 
the GB, with an energy of E^eg = —1.25 eV (c/. Fig. 2). 
A large part of this energy gain derives from ionic relax- 
ation during vacancy formation. Referring to Fig. 1, the 
high stability of vacancies in layer 1 is due to the close 
proximity of the pair of oxygen ions on either side of the 
interface and the electrostatic repulsion this causes. As 
one of the oxygens in the pair is removed, the other ion 
relaxes considerably, to a position in the symmetry plane. 

Having established the existence of low-energy vacancy 
positions close to the GB interface, we would now like 
to evaluate the corresponding equilibrium vacancy seg- 
regation and electrostatic potential barrier as function 
of temperature. To do so, we apply a space charge 
model"'^^"^^ with an ideal depletion approximation to two 
different cases: i) dry BZO, where a concentration of 
single-acceptor dopants are perfectly compensated for by 
double-donor vacancies, and ii) wet BZO, where disso- 
ciative absorption of water molecules from a humid at- 
mosphere is taken into account. 

In both cases we consider a one-dimensional model, 
with three regions of constant defect concentration: i) 
the GB core (0 < x < xq), where an increased concen- 
tration Cy of oxygen vacancies is expected and the pro- 
ton concentration is assumed to vanish, ii) a compensat- 
ing space charge layer (SCL) {xq < x < x*), depleted of 
oxygen vacancies and protons, and finally iii) the neu- 
tral grain interior {x > x*) with vacancy concentration 



Cy and in the wet case proton concentration Cqjj. The 
dopant concentration ca is taken to be fixed and equal 
throughout all regions and a concentration correspond- 
ing to occupation of 10% of the Zr sites is used in all 
calculations. Motivated by the high stability of oxygen 
vacancies in layer 1 (c/. Fig. 2), we define a core region 
enclosing that layer and denote the segregation energy for 
those vacancies with E^^^. Vacancies in layers 2,3,4,. . . 
are treated as having zero segregation energy and the 
symmetry plane, layer 0, is assumed void of vacancies. 
For the core half-width we use the value xq = 2 A. With 
this description the density of active oxygen vacancy sites 
in the core becomes about half that of the grain interior: 
^ 0.579 iV', where N' = S/ag = 0.039 A'^. 
We begin with the dry case and consider electrochemi- 
cal equilibrium between oxygen vacancies in the GB core 
and grain interior for an ideal solution^^: 
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The electrostatic potential difference between core and 
interior, A^'^ = (j){xo) — 4'{x*), is obtained by solv- 
ing Poisson's equation in the SCL with the boundary 
conditions 0(x*) = 4>'{x*) = 0. For the dielectric con- 
stant we use the value e = 32 eo^^ and the size of the 
SCL is determined from charge compensation with re- 
spect to the core: (2ecy — cca) xq = ecA {x* — xo), or 
X* — {2cy/ca) xq. It follows that the present model is 
applicable for Cy > ca/2 only. The solution for the po- 
tential difference is: 
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Eqs. 1 and 2 can now be solved itcratively to obtain Cy 
and A^'^ as function of temperature. Using the neutrality 
condition, the grain interior oxygen vacancy concentra- 
tion is in the present, dry, case given by: Cy = ca/2. Our 
numerical results are shown in Fig. 3 as dashed lines for 
a few different values of the segregation energy E^^^. 

Next, the wet case. Hydration of the grain interior is 
modeled by the hydration reaction (Kroger- Vink nota- 
tion) : 



H20(g)+Vo+OS^2 0Hc 



(3) 



Experimental values from the literature'^ are assigned to 
the hydration enthalpy AH^^^^. = —0.82 eV and entropy 
A'^hydr = —0.92 mcV/K and a water partial pressure of 
PH2O = 0.025 bar is used in the calculations. The interior 
hydroxide concentration Cqjj is obtained from the law of 
mass action, charge neutrality and site restriction'^: 
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where k ^ pa^o K and 
K = cxp{AS^^^JkB)exp{- 

die panel of Fig. 3 we show the result for the interior 
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FIG. 3. Core vacancy site-concentration and poten- 
tial barrier for vacancy segregation energies i?aeg ~ 
{—0.25, —0.75, —1.25} eV, where the thick hnes correspond 
to —1.25 eV. Dashed/sohd hnes have been used for the 
dry/hydrated case and dotted hnes indicate the bound of the 
present space-charge model, i.e. Cy > ca/2. The middle panel 
shows interior vacancy- and hydroxide concentrations for the 
hydrated case. 

hydroxide concentration Cqhi indicating the transition 
between wet and dry grain interior around 900 K. 

We are now in the position to determine the concentra- 
tion of core vacancies Cy and the electrostatic barrier Acf)"^ 
as function of temperature in hydrated samples. In this 
case Cy = {ca — Cqjj)/2 and TV' — Cy in Eq. 1 is replaced 
by TV'- 

'^OH ^ ^v- 1 ^iid 2 are solved together with 
the expression for the hydroxide concentration in Eq. 4 
and our final results for Cy and A(f)'^ in hydrated samples 
are shown as solid lines in Fig. 3. 

It is clear from Fig. 3 that, according to our model, 
segregation energies on the order that we find in the S3 
tilt GB can lead to significant electrostatic potential bar- 
riers as a result of oxygen vacancy aggregation at the 
GB. In dry BZO at 600 K, 21% of the oxygen sites hav- 
ing a vacancy segregation energy of —1.25 eV in the GB 
core are predicted to be vacant, giving rise to an electro- 
static potential barrier of 0.55 V. In the hydrated case, 
vacancies persist at the GB to the degree that, although 
the interior vacancies have been annihilated at 600 K, 
19% of the core sites are still vacant and the barrier is 
0.45 V. Even if the sample is cooled to room temperature, 
the number (assuming thermodynamic equilibrium) only 
drops to 17% and a potential barrier of 0.33 V is still 
present. 

In recent work, Kj0lseth et al^ used a space-charge 
model to evaluate the potential barrier corresponding to 
experimental results on grain interior and grain bound- 
ary proton conductivity in Y-doped BZO^'^"*'^^. A bar- 
rier on the order 0.4-0.6 V was found to match measure- 
ments made at 200-300°C. This agrees relatively well 
with our results, where a barrier of 0.4-0.45 V is seen for 
-^seg = -1.25 eV in this temperature range. Similar po- 
tential barriers fitted to experimental conductivity data 
were also found by Chen et al^" for 10 mol% Y-dopcd 



samples^'^'"'^°'^^ and by Iguchi et al^. Interesting to point 
out is the significant difference in potential barrier cal- 
culated from conductivity measurements of Duval^^, of 
BZO annealed at high temperature and of convention- 
ally sintered samples. Chen et al calculated a barrier of 
0.77 V in the conventionally sintered material, while the 
calculations of Kj0lseth et al show a barrier of 0.46 V in 
the annealed sample. 

When comparing our results to barriers estimated from 
experimental data it is important to keep in mind that 
our results are based on DFT calculation on one single 
grain boundary structure, while experimental results cor- 
respond to an average effect in polycrystalline samples. 
Theoretical investigations of additional GBs are there- 
fore warranted to conclude whether the magnitude of the 
present oxygen vacancy segregation energy is a general 
trait of GBs in BZO. Furthermore an aggregation of ac- 
ceptor dopants at GBs has been observed in BZO^'®'^®, 
an effect which we do not take into account and which 
would tend to diminish the positive charge due to oxygen 
vacancies. In addition, there is reason to investigate the 
segregation energy of protons in GBs, particularly con- 
sidering that the stability of protons in solid oxides has 
been seen to increase on oxygen sites where vacancies are 
also increasingly stable^'''^^. 

In conclusion, we present DFT calculations of oxy- 
gen vacancy segregation to a S3 tilt grain boundary in 
BZO. Using a space charge model we demonstrate that 
at 600 K, where the sample is fully hydrated in a wet 
atmosphere, the oxygen vacancies in the grain boundary 
core are still present with a site concentration of about 
20%. This corresponds to a Schottky barrier height 
of 0.45 V and compares well with recent experimental 
data®'®'^°'^^'^^. Our results indicate that oxygen vacancy 
segregation can be responsible for the low grain boundary 
proton conductivity in BZO reported in the literature. 

Since this paper was submitted for publication further 
relevant work on space charge and blocking effects at 
grain boundaries of BZO has been published^^^'^^ . 
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